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Abstract 

Following J.D. Murray, we consider a system of two differential equations that models traveling 
fronts in the Noyes-Field theory of the Belousov-Zhabotinsky (BZ) chemical reaction. We are 
also interested in the situation when the system incorporates a delay h > 0. As we show, the BZ 
system has a dual character: it is monostable when its key parameter re (0, 1] and it is bistable 
when r > 1. For h — 0, r + 1, and for each admissible wave speed, we prove the uniqueness of 
monotone wavefronts. Next, a concept of regular super-solutions is introduced as a main tool for 
generating new comparison solutions for the BZ system. This allows to improve all previously 
known upper estimations for the minimal speed of propagation in the BZ system, independently 
whether it is monostable, bistable, delayed or not. Special attention is given to the critical case 
r — 1 which to some extent resembles to the Zeldovich equation. 

Keywords: Belousov-Zhabotinsky reaction; comparison solutions; minimal speed; sliding 
solution method; bistable; monostable. 
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1. Introduction and main results 



One of useful objects associated with the famous Belousov-Zhabotinsky chemical reaction is the 



following dimensionless non-linear system 112 lL 122 1 



u t (t, x) — Au(t, x) + u(t, x)(l — u(t, x) — rv(t, x)), 

v,(t, x) = Av(t, x) - bu(t, x)v(t, x), ' 

called the Belousov-Zhabotinsky (BZ for short) reaction-diffusion system. The coefficients r, b 
are positive and u, v correspond to the bromous acid and bromide ion concentrations respectively. 
The front solution (u, v) = (0, 6){v ■ x + ct) of system ([TJ provides an appropriate mathematical 
tool for the description of planar waves propagating in a thin layer of reactant solution filled 
in a Petri dish [22]. Due to the chemical interpretation of ([TJ, only non-negative fronts are 
meaningful. Another requirement is the existence of the limits (0, S)(—o6) = (0, a), (cp, ff)(+oo) = 
(1,0) with a > 0. The exact value of a is not relevant: after rescaling u, v, we can take a = 1. 
By the experimental data 112 ll 12211 . r > 1. Nevertheless, almost all previous analytical studies 
of wavefronts (with two exceptions given in Propositions [TJ ITOb considered the case re (0^ 
which was proved to be of the monostable type. We observe that the standard definition 
of monostability/bistability needs an obvious modification in order to be applied to system ([TJ 
which has a continuum of non-negative equilibria. The degeneracy of the equilibrium (0, 1) is 
a special feature of model ([TJ complicating its analysis. For example, the recent Liang-Zhao 
general theory [ 14] of spreading speeds for abstract monostable evolution systems can not be 
employed here despite the fact that system ([1) is formally monostable and monotone for r < 1. 
This obliged us in B24ll to present a complete proof of the existence of the minimal speed of front 
propagation in ([TJ when r < 1. On the other hand, we show here that, for each r > 1, the BZ 
system possesses a unique wavefront solution, in full accordance with its formal bistability. 

Now, as it was argued in l24ll . a better theoretical prediction for propagation speeds in model 
([TJ can be also obtained by taking into account delayed effects during the generation of the 
bromous acid. For simplicity, and in order to connect with various analytical investigations, we 



will use here the following delayed version of ([T) proposed by Wu and Zou in [27]: 

u t (t, x) — Au(t, x) + u(t, x)(l — u(t, x) — rv(t — h, x)), 
v,(t, x) = Av(f, x) - bu{t, x)v(t, x). 

During the last decades considerable efforts have been made in studying the wave propaga- 
tion in ([TJ, (f2]) . The attention was focused on the stability, numerical approximation lfl9l I22I l23ll 

and existence e e h m m m m s m m of fronts. After linear changes, systems 



([TJ, <EJ acquire good monotonicity properties: they are quasi-monotone as partial differential 
equations H20L 12611 and they are monotone in the sense of Wu and Zou B27I1 . Hence, the front 
existence may be handled by the standard comparison technique well established for several 
decades 1I26I [27ll . Thus the existence of fronts for the BZ system is not longer an issue, in differ- 
ence with the determination or satisfactory approximation of the minimal speed of propagation 
in ([TJ, (f2j. Precisely this problem is our main concern here. It is quite noteworthy that a similar 
question (formulated as linear versus non-linear determinacy of the minimal speed) for a Lotka- 
Volterra reaction-diffusion competition model has received a considerable attention during the 
last few years . Finally, our secondary concern is the uniqueness of wavefronts (cf. (ll]): 

since these have to be monotone, we prove their uniqueness in the non-delayed non-degenerate 
case (i.e. r + 1, h = 0) by applying the Berestycki-Nirenberg sliding solution argument [|2|,[3|]. 
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1.1. Some previously known results 

For the sake of completeness, we state the most relevant known existence results for (H), ©. 
First of them was proved in 1I21II22I1 . it gives a lower bound for the admissible front speeds. Set 

c, := max [2% Vl - r, ( y/r 2 + 2b/3 - r)/ V2fc + 4r} . 



Proposition 1. Lef r,b > 0. If system (0 nas a positive componentwise monotone wavefront (or, 
shortly, monotone wavefront) connecting (0, 1) with (1,0) then c > q. 

It is easy to see that the estimation of Proposition[T]has the form c>c\-2 Vl - r forjrositive r < 
11/12 = 0.917 The next assertion summarizes the main existence results from IU2l[l3 , 281. 



Proposition 2. System (f7} has a positive monotone wavefront (u, v)(x, t) — (cf>, 8)(v ■ x + ct), 
\v\ = 1, connecting (0, 1) with (1,0) for each velocity 



c>c k = 



2 Vl - r, if rb + r< 1; 

2 V£, if either + r > l,b < 1, r € (0, 1] or = 1, r < 1; 
2, if > 1, rs (0,1]. 



Proof. The first condition was proved in [13, Theorem 3] under additional assumption r + b > 1 
when b < 1 . For r,b > satisfying rb + r < 1 it was also announced without proof as Theorem 



4.2 in [28]. The second and the third conditions were established in 1 13, Theorem 2]. □ 



Model (f2]i was considered in , from where we have the following 

Proposition 3. Assume that r e (0, 1). Then system (O has a positive monotone front (u, v)(x, t) — 
{<p,9)(v-x+ct), |v] = 1, connecting (0, 1) with (1,0) if either one of the following conditions holds: 
(I)b > 1 andc > max{b,2^}; (II) c > 2 Vl - r is such that b exp(-0.5c/z(c- yjc 2 - 4(1 ^ r))) + 
r < 1 . Finally, system (0 can nof nave wave/ronfi propagating at the velocity c < 2 Vl - r. 

Note that in the non-delayed case Proposition|3]is weaker than Propositions Q~|[2] 

Proof. See iflq. Theorem 3.1] for condition (/) and Ill7i Theorem 3.2] for condition (//). The 

final conclusion is known from [22] (for h — 0) and jl5ll (for h > 0). □ 



By 1121 1 Section 8], all wavefronts to ([T) are monotone. On the other hand, the delayed response 
may imply the loss of wave's monotonicity [JJl. Therefore it is worthy to emphasize that the 



inclusion of delay as in (|2]i does not change the monotone shape of fronts, see 11241 Theorem 6]: 



Proposition 4. If, for some r,b > 0, system (f2]) has a wavefront (u, v) — (<f>, 6) (v • x + ct), <p > 0, 
|v| = 1, connecting (0, 1) with (1,0), then <p'(t), -0'(t) > and 6(f), 4>(t) e (0, 1) for all t 6 R. 

1.2. General remarks about our approach and some useful relations 

The speed of front propagation in ([1), (f2]i can be estimated by means of the truncation method 
1 28], the shooting technique 012111311 and the upper and lower solutions I27I1 . Here, we 



use the latter approach complementing it by a useful idea about how to generate new comparison 
solutions. The main working tool will be regular super-solutions defined in Section l5TI Theorem 
[17] from the mentioned section is instrumental for the proofs of existence: its application with 
different regular super-solutions yields Theorems The same super-solutions are then used 
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in the bistable case, see Theorems [5] [9] Conceptually, Theorem [TTl is very close to highly non- 
trivial Theorem l(iv) from [3] (see also ||4j]). The proofs of TheoremfTTlancl the mentioned Chen 
and Guo result are, however, completely different. 

Asymptotic expansions of the eventual fronts at infinity are another key ingredient of our 
approach. In combination with a sliding solution argument they lead to 

Theorem 5. Let r + 1, b > 0. Then for each fixed admissible wave speed c the monotone 
wavefront (u, v) = (<p, 6)(v ■ x + ct), |v] = 1, connecting equilibria (0, 1) and (1,0) of system (|7J is 
unique (up to a translation). 

We also will need the following relations between the components of wavefront profile: 
Theorem 6. Consider <p, 6 as in Proposition^\and set ifr(f) := 1 — 6(f). We have 

A. Let r e (0, 1), K > \,L e (0, 1] satisfy K > b/(l - r) > L. Then 

I4(t - ch) < tff(t) < K(f>(t), t e R. (3) 

Ifb + r = 1, then <p(t - ch) < i[r(t) < tp(t). Hence, ifh = then <p = ij/. 

B. Let r>\. Then if/(t) > <p(f), t e M. 

C. Suppose that r e (0, 1], then i// 2 (t) < M(p(t), fel, M := max{l, 2b}. 

By part [A], the BZ system with h - 0, b + r - 1 essentially reduces to the KPP-Fisher equation 
J7l|2l|]. Part [B] has a clear chemical interpretation: the sum of the (normalized) concentrations 
of the bromous acid and bromide ion in the propagating wavefront is strictly less than the con- 
centration of the bromide ion far ahead of the wavefront. Part [C] connects (f2]) with the delayed 
Zeldovich equation u,(t,x) = Au(t, x) + f3u 2 (t - h,x)(l - u(t,x)), /3 = minjb, 0.5). Actually this 
relation suggested the correct form of asymptotic expansions (0 below (see also [24, Lemma 26 
and Corollary 27]). 



1.3. Main results: monostable case 

For the non-delayed BZ reaction (Q3 and r 6 (0, 1), the existence of the minimal speed of front 
propagation c*(Tl) was proved in 11261 p. 333]. The speed c*(Tl), however, is minimal only for the 
fronts taking values in special domains II called the balance polyhedrons. Since the BZ system 
has a continuum of equilibria, none of these domains can cover the whole region admissible for 
wavefronts, see [26, Fig. 5.1, p. 334]. The existence of the positive minimal speed independent 
on II was established in ll24l Theorem 7], by means of regular super-solutions. By Theorem 8 
below, c t - 2 Vl - r if rb exp(-2/z( 1 - rj) + r < 1. However, due to Proposition!]] it may happen 
that c t is not linearly determined (i.e. c* > 2 Vl — r), cf. ifsL I9I [loll . Even for the non-delayed 
BZ system, the exact value of c* in the case rb + r > 1 is unknown and represents an interesting 
open problem. The next theorems show that the use of regular super-solutions in the Wu and Zou 
approach [27] yields important improvements of the estimations of c* even for the non-delayed 
model. Set b' :- be~^ 2 and let c# = c#(r, b, h) be the unique positive root [24] of the equation 

c = 2maxkvT37,-^=Uf 2 ^ tf ^ V (-2h(l - r)) + r < 1; 

1 VT+17| I 2VF/VT+17, if rbexp(-2h(l-r)) + r > 1. 
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Theorem 7. Let r e (0,1], c > c#. Then system (0 has a positive monotone front (u,v) 
(4>, 6)(v ■ x + ct), \v\ — 1, connecting (0, 1) w/fft (1,0) and such that (i) if r = 1, c > c#, ffte/1 



2c 2 /ft 8c , 1 ln(-f) 1 

t L 3ft ft t s P 



0(t) = l + ^-^(c 2 (l + ft+i)-4)^^ + 0(i), r^-c»; 
f 3 b t t 

(ii) if r e (0, 1), c > 2 Vl — r, then, for some s > ant/ A :— 0.5(c — -^c 2 — 4(1 — r)), /f holds 

U e M}-ch) 

(p(t) = e At + 0(e {x+E)t ), 0(f) = 1 - — + 0(e {A+E)> ), t -oo. (6) 

1 — r 



Theorem 8. Assume that r e (0, 1], c e [2 Vl - r, c#) ant/ 

t (ii* h c 2 w, 1 — r 

/(c 2 , r, ft, A) := c 2 (-^ + -) + In — - > 0, (7) 

8r 2 4or 2 r 

where w, = 8.21093 . . . denotes the greatest positive root of the equation a> - A + 2\na>. Then 
system (0 has a positive monotone front connecting (0, 1) with (1,0). Asymptotic formulas (|3J 
(when c > 2 Vl — r) and ((5]) (when r — 1) are fully applicable for this wavefront. 

Observe that inequality (0 can be written as c > c D = c a (r, b, Ii) where c a is the unique positive 
root of the equation f(c 2 , r, b, h) — considered with fixed r, b, h. 

1.4. Main results: bistable case 

The next assertion can be considered as a dual to Theorems [7J [8] Indeed, it essentially 
amounts to the non-existence of bistable waves for c > c#(r, b, h) and c > c a (r, b, h): 

Theorem 9. Let r > 1, b > 0. Then system (0 has at most one (a unique, if h — 0) positive 
monotone wavefront (u, v) = (cp, 8) (v ■ x + c+i), (f> > 0, |v| — 1, connecting (0, 1) with (1,0). The 
( unique) velocity of propagation c* satisfies the inequality c*(r, ft, h) < min{c#(r, ft, h), c a (r, ft, h)}. 
In addition, c*(r, ft, h) is non-increasing in h. 

The wave existence problem for the bistable BZ delayed system requires a different approach 
and it is not considered here. In the non-delayed case, the wavefront existence was established 
by Kanel in [13, Theorem 4]. In view of Theorem[9] Kanel's result can be reformulated as 

Proposition 10. Let h — 0,r > 1. Then system (0 has a positive monotone wavefront for the 
speed c* such that c K :- b/(2 V(r + ft) [min( 1 , b)(r + ft) - 0.5ft]) < c+ < 2 Vmin(l,ft). 



Table 1 : Analytical and numerical estimations of c» , 





Propositions |2J 1 1 1 


Theorems|7J|2l 


Theorems|8]|2] 


Numerical c. 


Propositions LjJI 1 0| 


(0.5; 5) 
(0.5; 10) 

(i;5) 
(5; 0.5) 


c> 2 
c> 2 
c>2 
c* < 1.41... 


c> 1.82... 
c> 1.90... 
c> 1.82... 
c* < c # = 1.15... 


c> 1.62... 
c> 1.71... 
c> 1.47... 
c* < Co = 0.59 . . . 


c, rj 1.46... 
c» rj 1.50... 
C RJ 1.13... 

c*»0.12... 


c. > ci = 1.414... 
c. > ci = 1.414... 
c, > a = 0.289 . . . 
c K = 0.067..., c, = 0.007... 
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Example In Table 1, for h — 0, we compare results of Theorems |7]|8]|9] with previously known 
ones (Propositions l2l [Toll. Notation like c > 1.82 .. . means that system ([TJ has a positive front 

for each velocity c > 1 .82 Numerical estimations of the minimal speed c* are taken from 

11191 Table 3] and [23, Table 1]. Lower bounds for c», c* are computed from Propositions [T1 [TOl 

Finally, the organization of the paper is as follows. Sections l2l l4l 15 .21 15 . 31 and l6l contain the 
proofs of Theorem [6] [5] Q [8] and [9] respectively. Asymptotic behavior of profiles at infinity is 
analyzed in Section[3] Our main technical result (Theorem [TTb is proved in Section lBTTI 



2. Proof of TheoremE] 

Let (m, v) = ((p, ff)(v ■ x + ct) be a wavefront to (0. After introducing 0(f) = 1 — 0(t — ch), we 
obtain the following boundary value problem for the determination of fronts in the BZ system: 

0"(f) - ccf>'(t) + 0(/)(l - r - 0(f) + nftt)) = 0, 

0"(f) - c0'(f) + b<f>(t - ch){\ - 0(f)) = 0, (8) 

> 0, < 1, 4>(-oo) - 0(-oo) = 0, 0(+oo) = t/r(+oo) = 1. 

[A] Set z(t) := K<p(i) - if/(t). It is easy to see that 

z"(t) - cz'(t) + {K(l - r)(f>(i) - b<f>(t - ch) + K<t>(t)(n/,(t) - (p{t)) + bif/(t)<p(t - ch)} = 0. 

Since z(— 00) = 0, z(+oo) = K — 1 > 0, the non-positivity of z at some points implies the 
existence of some r such that z(t) < 0, z'(r) = 0, z"(t) > 0. But z(t) < implies K(f>(r) < 0(t) 
and therefore 

= z"(t) + {K(l- r)4>(r) + K4>(T)(rif/(T) - 0(t)) + b(i//( T ) - l)0(r - ch)} 

> {K(l- r)<p{r) + K<p(r)(np(r) - 0(r)) + 2#(r) - 1)0(t)} 

> 4>(t) {[K(l -r)-b]+ Kcf)(T)(rK - 1 + b)} > 0, 

a contradiction. The latter inequality holds obviously if rST — 1 + b > 0. If rK — 1 + b < 0, then 

#(1 - r) - b + K(f>(T)(rK - 1 + b) > K(l - r) - b + ZT(r/: - 1 + b) = (rK + b)(K - 1) > 0. 

Next, set z(t) := L$(t - ch) - t//(t). We have z(-oo) = 0, z(+oo) = L - 1 < 0, so that the non- 
negativity of z at some points would imply the existence of some r such that z(j) > 0, z'(j) = 0, 
z"(t) < 0. But then L(/>(t - ch) > 0(t) and therefore 

= z "(t) + (L( l-r)- b)4>(T - ch) + Lc/)(t - ch)(rif/(T - ch) - 4>(t - ch)) + b^(r - ch)t//(T) 

< (L( l-r)- b)(p(T -ch) + L(p(T - ch)(riJ/(T) - <p(r - ch)) + b<p(r - ch)i//(r) 

< L, := (L(\ - r) - b)(f>(T - ch) + L0 2 (t - ch)(rL - 1 + b) < 0, 

a contradiction. The latter inequality is obvious if rL - 1 + b < 0. If rL - 1 + b > then 

U < <p(r - ch)((L(l -r)-b) + L(rL - 1 + b)) = 0(r - ch)(rL + b)(L - 1) < 0. 

[B] Consider z(t) := 0(f) - (f>(t). We have that z(+oo) = 0, 

z"(t) - cz'(t) + b(p(t - ch)(l - 0(f)) - 0(0(1 ~r- 0(f) + r0(f)) =0, f e E. (9) 
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If z(s) < at some s then there exists r such that > z(t) = min, E R z(t). We have that z"(r) > 0, 
z'(t) = 0, <Mt) < 0(t), and 

-0(t)(1 - r - 0(t) + njf(r)) > -0(t)(1 - r - 0(t) + nflr)) = 0( T )(r - 1)(1 - 0(r)) > 0, 

contradicting to ||9). 

[C] Consider z(f) := M<f>(t) - if/ 2 (t). Since M > 1, we have z(-oo) = 0, z(+oo) = M - 1 > 0. 
Thus the non-positivity of z implies that z(t) < 0, z'(r) = 0, z"(t) > for some rel. Hence, 

M<P(t) < </r 2 (r), M0'(t) = 2^(t)^'(t), M0"(t)>2(A(t)i/t"(t) + 2(i/''(t)) 2 , 

> 2if/(T)f(T) + 2(i//'(t)) 2 - 2a/r(T)f(T) + M0(t)(1 - r + r^r) - 0(t)), 
= 2^(t)(/t"(t) - 2o//(t)i//'(t) + 2^(t)0(t - c/i)(l - ^(t)), 

so that 

> 2(\p'{T)f + M4>(t)(\ - r + np(T) - 4>(t)) - 2b\l/(T)4>(T - ch){\ - iff(r)) 

> M<p( T )(\-r + r<A(r) - 0(r)) - 2^(t)0(t)( 1 - -A(t)) 

> 0(r) {M(l - r) + tff(r)(Mr - 2b) + t// 2 (T)(2b - 1)} > 0, 

a contradiction. Here we observe that the polynomial p(z) := M(l —r) + z(Mr- 2b) + z 2 (2b - 1), 
z := >J/(t) € (0, 1), satisfies p(0) = M(l - r) > 0, = M - 1 > 0, so that p(^(t)) > if 
2b - 1 < 0. If 2b - 1 > then we choose M := 2b > 1 to obtain 

p(z) = 26(1 - r) - 2fe(l - r) + z 2 (2Z? - 1) = 26(1 - r)(l - z) + z 2 (2Z? - 1) > 0. 



3. Asymptotics of wavefront profiles 

First, we observe that the derivatives (/>', ift' of wavefront components are bounded and uni- 
formly continuous on R so that 0'(±oo) = i/r'(+oo) = 0. This fact is well known (cf. [27, Section 
2]) and its proof is omitted. Incidentally, the relation i[f'(±oo) = implies the positivity of each 
admissible speed (i.e. c > 0): it suffices to integrate the second equation of © on R . 

Next, assume that r e (0, 1]. Using Theorem[6lB] if r — 1 and integrating ([8]l on (-oo, t], t < 
t a , we get, for sufficiently large negative t a , 



<f>'(f) <<f>'(f)+ f <t>(s)(l-r+ri}r(s)-<f>(s))ds = c(f>(t), <y(t)+b f (f>(s-ch)(l-i/f(s))ds = ctff(t), 

and therefore z(f) :- <p'(t)/<p(t) < c, t < t a , <p € Lj(R_). Furthermore, z satisfies the equation 

z' + z 2 - cz + (1 - r) = /(f), where /(f) := 0(f) - n/Kf). (10) 
Let A = /1(c) < ju = /i(c) denote the roots of the characteristic equation x 2 - cx + (1 - r) — 0. 



Lemma 11. Let (cp, ip) be a traveling front q/dS]) and r e (0, 1]. Then (a) c >2 Vl - r, (b) there 
exists finite limit lim (f>'(f)/(f>(f) e as t — > -oo. 

Proof. Recall that /(f) — > as f — > -oo. (a) Suppose that c < 2 Vl - r. Then, for some tb < t a , 
it holds /(f) - z 2 + cz - (1 - r) < for (f, z) € (-oo, f fe ] x [0, 2c]. However, as a simple analysis 
of the direction field for equation ([TUt shows, this contradicts to the property z(f) e (0, c), t <t a . 
(bl) Let c = 2 Vl - r and take some small e > 0. By analyzing the direction field again, we 
can see that there exists t c such that (f, z(f)) e (— °o, ? c ] x ( _e + c/2, c/2 + e) for f < f c . Hence, 
z(f) — > c/2 = /I = yu as f — > -oo. (b2) The situation when c > 2 Vl - r is similar to (bl). □ 
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Corollary 12. Let r e (0, 1). Then there are t\, m e {0, 1} and v(c) e {A(c),p(c)} such that 
(tffl + fi), <p(t + h), (j)'{t + h)) = (-t) m e v(c) '(be- v(c)ch /(l - r), 1, v(c))(l + o(l)), t -> -oo. 



Proof. By Lemma fTTl (f>(t),(p'(t) decay exponentially at -oo. Then <A(f) has the same property 



due to Theorem[6]A]. Therefore we can apply Proposition 7.2 from 1 18] together with Theorem 



[6[A] to system ([8]) in order to obtain the above asymptotic formulas for <p, <p', ift. Note that m = 1 
only when c — 2 Vl - r. □ 

Lemma 13. Let (<f>, if/) be a wavefront for and r > 1. 77zen, /or some A > 0, ?2 € R, a/id 
sma/7. cr > 0, /f ZioZcfc (f>(t + t 2 ) = e^ c)t + 0(e ( - 2c -< f)t ), if/{t + t 2 ) = Ae c > + 0(eW c >- <r >'). t -co. 

Proof. Integrating the first equation of © from -co to t, and using the inequality 1 — r — (p(t) + 
rcp(t) < for all large negative t (say, for t < T where, simplifying, we can take T = 0), we obtain 
that 0'(O _ c< X f ) > for f < 0. Thus (f>(t) < (f>(0)e cl , t < 0. Similarly, from the second equation of 
(O, we deduce if/(t) > if/(Q)e cr , t < 0. The latter equation can be written as if/"(t) - cif/'{t) = F(t), 
where F(t) := b<f>(t - ch)(if/(t) - 1) = 0{e ct ), if/(t), if/{t) = o(l), t -co. But then HI Proposition 
7.1] guarantees that iff (t), iff' (t) = 0(e (c ~°~^'), t — > -co, for each small cr > 0. Now, writing the 
first equation of © as <j>"{t) - c(f>'{t) + (1 - r)cf>(t) = G(f), where G(f) = CKe* 2 ^'), f -» -co, 
we find analogously that <f){i) = Be M(c)r + 0(e (2c_cr) '), t -> -co, where cr > 0, B > 0. To prove 
that B > 0, it suffices to repeat the proof of Lemma QT] (note that z(t) is bounded on R because 
otherwise it blows up in a finite time). Hence, F(t) - 0{e (M(c ^), t — > -co, if/(t) > if/(Q)e cl , t < 0, 
p(c) > c. By 118, Proposition 7.1], this yields the required asymptotic formula for if/. □ 

Next, we consider the case when t — > +co. In order to linearize system ((8} along the positive 
steady state (1, 1), we use the change of variables (pit) = l-^(f), if/(t) = 1 - 0(t - ch), which leads 
to 

(0 - (f) - #0(1 - f (f) + r0(f - c/z)) + r0(f - cA) = 0, 
6>"(f) - c6»'(f) - b6(t)(l - £(?)) = 0. { ' 

The characteristic equation (z 2 — cz - l)(z 2 - cz — b) = for this system at the zero equilibrium 
has two positive (£2, £2 = 0.5(c + Vc 2 + 4fe)) and two negative eigenvalues (£ 1 and (\ = 0.5(c - 
Vc 2 + 4b), respectively). 

Lemma 14. Let r > 0. Then for some appropriate A > 0, ?o, c/, rfi ana? sma// cr > 0, we /iave f/iaf 
(t/>(t+t ),<t>'(t + t )) = -Ae^(l,?i)+ 

f (1 - reM'- ch) /(b - 1), -r^eM'-^Hb - 1)) + 0(e^-°">), & # 1, 

1(1- r(f + t/)ef' ( '- c/,) /(c - 2ft), -rft(* + di)ef' (f - c/l) /(c - 2ft)) + Oie^-^'), b = I, 



(if/(t + t Q ), if/(t + t )) = (1 - e M '- ch> , -^e M '- ch) ) + Oie^-^l t 



+ CO. 



Proof. Since #(+00) = £(+co) = and the linear system /'(f) - cy'(t) - y(t) + rz(t - ch) = 0, 
z"(t) - cz'(t) - bz(t) = possesses an exponentially dichotomy on R + , the perturbed system 

y"(t) - cy'(t) - y(t)(l - ((t) + r6(t - ch)) + rz(t - ch) = 0, z"(t) - cz'(t) - bz(t)(l - £(t)) = 

is also exponentially dichotomic on R + . As a consequence, we obtain that 8{t),ff (t),^(t),^'(t) = 
0(e 1 '), t — > +00, for some negative I. Moreover, by applying the Levinson asymptotic integration 
theorem 16|] to the second equation of (TTTb , we find (cf. [7, Lemma 19]) that, for some to, 

(6(t + to), &{t + to)) = (e^"(l + o(l)), -fte^U + o(l))), t -> +co. 



Then 018 



Proposition 7.2] applied to the second equation of (fTTb yields the required estimation 



(G(t + t ), 8'(t + f )) = (e^', fte fl< ) + 0(e® _cr)f ), t -> +00. (12) 

Let simplify O by assuming t = 0. If b + 1 then ft * ft andy = ^(t)+re il{ '- ch) /(b- 1) = 0(e") 
satisfies 

y"(0 - cy'(0 - y(t)(l + m(t)) = 0(e ill -°' )t ), t -> +00, (13) 



where m(f) = 0(e ). Applying again Proposition 7.2 from II 1811 . we conclude that if ft > ft 
(equivalently, o e (0, 1)) then y(t),y'(t) = Oie^''""") with cr' e (0,cr) so that 

(m?(t)) = + 0( e »-'>'), f -» +00. 

l-o 

When ft < ft (that is b > 1), we find similarly that, for some A > and t — > +00, 

< f(f) = Ae ?1 ' + r^ l( '- cA) /(l-fe) + 0(c (fl -' r ' )f ), £'(0 = Afte ?1 ' + rft/ l( '- c ' !) /(l-») + 0(e &_!r ' ) ')- 
Finally, if = 1 then ft = ft and therefore 

2ft -c 

satisfies (1131 . As a consequence, we obtain (once more invoking Proposition 7.2]) that, for 
some real d, d 2 , it holds (ftf)(0 = r(f + rf, ftr + d 2 )e M '- ch) /(c - 2ft) + 0(e (f '- </) '), f -» +00. □ 



4. Proof of Theorem|5] 

The proof is based on the Berestycki-Nirenberg sliding solution argument. Let (0i, 00, (02, 02) 
be two different (modulo translation) traveling fronts of ((8} considered with h — 0. By Lemma 
[l~4l without restricting the generality, we may assume that 0i and 02 have the same first terms of 
their asymptotic expansions at +00. In addition, due to Lemma [T3l (employed when r > 1) and 
Corollary[T2](for r < 1), we can index iftj in such a way that either 0i(f) > 02(0 on some infinite 
interval (-00, T] or 0i, 02 also have the same first asymptotic exponential terms at -00 (recall 
that r + 1). In each case, the closed set >S := \s : if/\(t + s) > 02(O> t e R] + R is non-empty and 
contains finite s* :=inf»S. Similarly, there exists the leftmost t m such that <p\(t+t t )> fait), tel. 

Let us show that actually s* - 0. Indeed, if s* > then, due to the chosen asymptotic 
behavior of ij/j at ±00, we find that, for each e e [0, s»), it holds 0i(f + s* - s) > $2(1) for all 
I e I excepting t from some compact interval. This implies the existence of finite t such that 
8(f) = 0, 5"(t) > 0, 5(t) := ifri(t + s t ) - 02(0 ^ 0. If we suppose additionally that s» > f* 
then 0i(f + 5*) - 02(0 > 0, f G R. (Note that 0i(F + s») - 02(0 = implies that s» = ?* and 
<p\(t + s t ) - <p' 2 (t) = 0. Since also 5(t) — = d>'(f) = 0, the solution uniqueness theorem for © 
assures that (0i , 0i)(? + s») = (02, 02)(O)- But then we get from (O the following contradiction: 

= 5"(t) - c6'(t) + b(Mt+ *.) - 02©)(1 - 02(0) > 0. (14) 

Hence, we have to consider the case when t m > s* > and 0i(f + s*) - 02(0 ^ 0. Note that 
6(±oo) = 0,(5(0 > 0, and therefore 6(t) has at least two local maxima at some tf. t\ <t<t%. 
Since 6"(tj) < 0, S'(tj) = 0, estimations similar to (fL4l > shows that (f>i(tj+s»)—<fa,(tj) >0,j- 1,2. 
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Next, set S a (t) := fa(t + s* + a) - fa(t). Functions S a (t) are increasing in a and strictly positive 
on [fi,?2] for all large a > 0. On the other hand, So(t) has at least one zero on (fi,?2)- This 
means that for some a t > and t c e (t\, t<i) function S*(f) :- S a Xt) reaches at t c its zero global 
minimum on [t\, f 2 ]. Therefore 5" (f c ) > 0, 5»(f c ) = 0,S„(t c ) = 0, so that, due to ©, 

= S'J{t c ) - cS'Jt c ) + rfa(t c )(fa(t c + s t + a,) - if/ 2 (t c )) > 0. (15) 

This shows that a* — and that 

fa(t + s t ) - fa 2 (t) = fa(t + s.) - ^ 2 © = + s.) - 02© = 0i( f ~ + s *) - 02© = 0. 

But then, by the uniqueness theorem for ©, (fa, fa)(t + s*) = (fa, fa)(t), t e M contradicting to 
our choice of [fa, if/j). Therefore we conclude that s* — and 5{t) > 0, t e M. In the remainder 
of the proof we will analyze three possible mutual positions of f* and 0. 

Case A : f, < 0. Recall that i^i(f) > (f^CO. + f ») ^ 02(0- Due to the coincidence of the 
principal terms of asymptotic representations for fa, fa at +oo, we see that, for every small 6 e 
(0, |?*|) the graphs of functions fa(t - 6) and fa(f) have at least one intersection on some interval 
[T, +oo). In fact, we may assume that fa(T - 6) > fa(T) and fa(t - 6) < fa(t),t € [T\, +oo), for 
some T\ > T. It is clear also that fa(t - 6) > fa(t) for all t e R. Next, we consider the family of 
functions fa (t - 5) + a and the following non-empty and closed set 

21 := {a > : fa(t - 5) + a > fa(f), t e [T, +oo)}. 

Set a* = inf 21, it is evident that a„ > and that w(t) := fa(t — 8) +a m — 4*i(t) has at least one zero 
t p e (T, +oo), where, in addition, w'(t p ) — 0, w"(t p ) > 0. But then, due to equations ([8]), 

= w"(t p ) - cw'(t p ) + b{a t fa(t p ) + [fa(t p -6)- fa(t p )](l - fa(t p - 6))) > 0, 

a contradiction proving that t m > 0. In fact, we have established a stronger result: for every 6 > 0, 
the inequality <p\(t - 5) > fa(t) does not hold on any infinite interval [T, +oo). As a consequence, 
there exists a minimal p e [0, f,] such that fa(t + p) > fa(t) for all t e [T, +oo). That is, for 
every small 6 > 0, equation 0i(f + p - 6) = fait) has at least one root on (T, +oo) (otherwise, 
fa(t + p - 6j) < fa(f), t > T, for some 5j — > and therefore 0i(f + p) < 02© t>T, implying a 
contradiction: fa(t + p) = fa(t), fa(t + p) > fa(t), t > T). 

Case B : t t = 0, so that fa(t) > fa(t), fa(t) > fa(t), t e R, and, for each 5k > 0, the inequalities 
fa(t - 5\) > fa(t), fa(t - 62) > fa(t) do not hold on any interval [T, +00). Now, it is easy to see 
that, in fact, S*(t) := fa(t) - fait) > 0, t e R. Indeed, otherwise S t (t c ) = for some t c and thus 
we get a contradiction as in (fTBT l. where s* — a* — should be taken. Hence, for a fixed T and 
for small 6 > 0, each difference fa(t — S) — fa(f), fa(t-5) — fa(t) has at least one zero on [T, +00). 
We can choose large T and small 6 > in such a way that 

fa(T)-2r(l-fa(T-5))>0, fa(T - 6) > fa(T), fa(T - 5) > fa(T) > 2/3. (16) 

In the next stage of the proof, we apply the sliding solution argument to the families e + fa(t - 5) 
and 2er + fa(t — 6). It is clear that the sets 

£1 := {e > : e + fa(t -5)> fa(t), t e [T, +00)}, 



£2 := {£ > : 2er + fa(t -5)> fa(t), t £ [T, +00)} 
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are closed and non-empty, and that ej = inf£, are positive. Suppose first that e\ > e 2 . The 
difference y(t) := e\ + ifti(t - 8) - if/ 2 (t) reaches its global minimum at some point t m > T where 
y(t m ) = Y(t m ) = and y"(t m ) > 0. We also have that 

2e\r + <p\(t m -S)> 2e 2 r + (p\{t m - 8) > (f> 2 (t m ). 

Therefore, using ((8} again, we find that 

= y"(t m ) - cy'(t m ) + b [ei4> 2 (t m ) + (4>i(t m -8)- 2 (U)(1 - #i(t m - 8))] > 

be x [fefcO - 2r(l - fr(t m - 8))] > be x [fa(T) - 2r(l - ifr^T - 8))] > 0, 

a contradiction. So e\ < e 2 and the difference a(f) :- 2e 2 r + <f>\(t — 8) - 02(0 reaches its global 
minimum at some point t„ > T where a(t„) = a'(t n ) = and a"(t n ) > 0. We also have that 

e 2 + ipi(t n -5)>ei + \{/i{t n -S)> iff 2 (t„). 

But then, after invoking (fSJ, we get a contradiction: 

= a"(t„) - ca'(t„) + rfaitn - 5)[2e 2 + \f/i{t n -8)- if/ 2 (t n )] + 

2e 2 r(r - 1 + (f>i(t„ -8) + 2e 2 r - nf/ 2 (t„)) > 2e 2 r(r(l - iff 2 (t„)) - 1 + 1.5<pi(t„ -8) + 2e 2 f) > 0. 

Case C : t* > 0. For a fixed large T > 0, we consider <p\{t + p) where p e [0, ?*] was defined in 
the last lines of subsection 'Case A'. Then ij/\(t + p) > i// 2 (t), t € M, and, for each small 8 > 0, the 
equation <pi(t + p - 8) = <p 2 (t) has at least one root on {T, +oo). From this point we can follow the 
proof given in Case B (beginning from ([Toll). Actually, it will be literally the same proof if p = 0. 
If p > we have to replace, starting from ([16) , (f>\(t-8), 4>\{t-8) with 4>\(t + p-8), \jj\{t + p — 8). 
Note also that e x = 0, e 2 > if 8 e (0,p). 

5. Regular super-solutions and proof of Theorems |7j|8] 

Assume that r > and c 2 > 4(1 - r). Recall that A = A(c) < p — p(c) denote the real roots 
of the characteristic equation x(z, c) :- z 2 - cz + (1 - r) = 0. Fix some positive v e (A.,p). If 
r G (0, 1) then we define k as the maximal positive integer such that kA < v and (k + I) A > v. 
Obviously, if k > 1 then we havexijA, c) < for all j = 2, . . . ,k. 

5.1. Regular super-solutions and a preparatory theorem 

To prove the existence of monostable fronts, we will use Wu and Zou version 112711 of the 
upper and lower solutions method. Below, we propose a trick which increases the effectiveness of 
this approach for the BZ system. We will show that it suffices to find only two solutions (instead 
of four ones which must agree amongst themselves) of a system of differential inequalities. 

Definition 15. Assume that continuous and piece-wise C 1 -smooth functions <p + are posi- 
tive and have positive derivatives in some neighborhoods 0\, 2 of the sets (-°°, fi], (— oo, t 2 ], 
respectively. We admit here that (if/' + , cf>' + ) has a finite set D — {d\ < d 2 < ... < d/^), d\j < 
min{fi,?2}, of the discontinuity points and one-sided derivatives of iff + , <p + satisfy ip' + {dj-) > 
>p' + (dj+), <p' + (dj-) > 4>' + (dj+). Suppose also that t/r + (-oo) = + (-oo) = 0, ift+ih) - 4>+(h) — 1, 
and that ifr+, (p+ are C 2 — smooth in some vicinities oft\,t 2 and that 
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Dl. For a fixed positive v 6 (A, p.), m e {0, 1), and some positive constants C\, e, it holds 

if, + (t) = 0(te vt ), (0 + (f),0' + (O,<(O) = C 1 (-f)'" e w (l,v,y 2 )(l t ^ -ex,. 

D2. /ff < minffj, f 2 ), f £ 2), then 

Ai(0 + , tA+XO : = 0i'(O - e#.(f) + - r - 4> + (t) + f^+CO) < 0, 



A 2 (0+, := Kit) - cif/+{t) + bcf> + (t - ch){\ - iff + (t)) < 0. 
D3. If t\ < t 2 then 0"(f) - ccf>' + (t) + (f>+(t){l - cf)+(t)) < 0, t e [t u t 2 ]. 
D4. Ifti > t 2 then if/"(t) - cif/'Jt) + />min{l, (f>+(t - ch)}(l - iff+(t)) < 0, t e [t 2 , h]. 



(17) 



We will call such {iff + , + ) a regular super-solution for (0. Observe that we may suppose that 
<f> + is defined, strictly increasing and smooth on [t 2 , +oo), this fact is implicitly used in D4. 

Remark 16. Suppose that <f>+,if/+ are increasing and that inequalities ( 1771 ) hold for all t < 
max{?i,?2). Then conditions D3,D4 are satisfied automatically. Indeed, in case D3, we have 
that Ai(0 + , 1) < Ai(0 + , iff+) < 0, t e [t\,t 2 ], while, in case D4, 

iff'lit) - ciff' + {f) + bmm{\,(f> + {t - ch)}(\ -if/+(t)) < A 2 (0 + ,iA + ) < 0, t e [t 2 ,h]. 



Note that the upper solutions for the BZ system proposed in earlier works (e.g. see 11171 12711 ') 
have 'correct' behavior at -oo and therefore do not satisfy condition Dl. 'Correct' here means 
'asymptotically similar to the true wavefront' (i.e. satisfying Q, ©). 

Theorem 17. Suppose that for given parameters b,c > 2 Vl - r, r 6 (0, 1], h > 0, system (0 
has a regular super-solution (iff + , (f> + ). Then there exists a monotone wavefront for (0 moving at 
the velocity c and satisfying (0, ((6)). 

To prove Theorem [T7] we will need several auxiliary statements. The first of them can be 
viewed as a variant of the Perron theorem for piece-wise continuous solutions, cf. flO]. 
Lemma 18. Let iff : R — > R be a bounded classical solution of the impulsive equation 



iff" + Aif/ + Biff - /(f), A#, = a h Aiy\ tj = B h 



(18) 



where {tj) is a finite increasing sequence, f : R — » R is bounded and continuous at every t + tj 
and Aw\ tj := w(tj+) — w(tj—). Assume that ^\ < < £, 2 are real roots ofz 2 + Az + B — 0. Then 



— — if e^'- s) f(s)ds+ f e^'- s) f(s)ds 

1 - ?2 \ J-oc Jt 

X*'*'"" -ft) + J> fl( '~ 0) ^ a > -ft) 



(19) 



Proof. It is straightforward to check that iff defined by ( fT9l verifies equation (1181 1. 
Lemma 19. For r e (0, 1), set a\ := 1, b\ := be~ Ach f(\ - r). There are functions 

and a polynomial P(x, y) such that, for all t e R, 

<f> A \t) - ccf>' A (t) + A (O(1 - r - M0 + nff A (t)) = A 2 P(e A ',A)e A(k+l) \ 
iff'^t) - cif/' A (t) + b(f> A (t - ch) = 0. 
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t + t: 



□ 



A(a x e M + a 2 e 2Ar + 
A(b x e At + b 2 e 2A ' + 



+ a k e kA >), 
+ b k e kA >), 



(20) 



Proof. Indeed, for a suitable polynomial P(x, y), we have that 
<(f) - c<p' A (t) + <p A (t)(l - r - <p A {t) + rifr A (t)) 

k 

= A YfxUX, c)aj - J] Aa p {a q ~ rb q ))e^' + A 2 P{e A ' , A)e A{k+l)t , and 

7=1 p+q=j J^-. 

^'(t) - af,' A (t) + b(/> A (t - ch) = A 2Jx(j*, c)bj + (ba je -i Ach + (r - l)bj)]e^. 

7=1 

In order to obtain (l20l >. we define recursively (J - 2, ... ,k) 

be- Ach A ~L P+q =ja p (a q -rb q ) ba f -i Mh 

a\ = l.b] — , a, — A , b i — . U 

l-r' J xij^c) 1 l-r- X (M,c) 

Remark 20. It is easy to see that, for some rational functions aj(b, r, A, c, t) and b jib, r, A, c, t), 
it holds 

aj = A'- X a :i {b, r, A, c, e - Ach ), bj = A hX bj{b, r, A, c, e^' 1 ). 

Therefore A-\(/) A (t), if/ A (t)) = {l,be- Ach (\ - r)- l )e A ' + AO(e 2At ), t -> -oo. This implies that for 
every to e R there exists Aq > such that (p A ,i// A ,<p' A , \j/ A are positive for all t < to, A e (0, Aq]. In 
addition, the derivatives of(p A , \j/ A have the property lirrj^o+wf (f), ifr. (t)) — (0, 0), k = 0, 1,2, 
uniformly on (— oo, to]. 

Next, in order to get an analog of <p A , i/r A when r = 1 , we consider the functions 

2c 2 lb yUn(-f) rin 2 (-f) 
Mt) : = __ + __ + __, 

2c Cln(-f) r Q\n 2 {-t) 
:= ~T + ^^ + 7 + ^^> l< - e > 

which coefficients C, F depend only on c, b, h and are defined explicitly by 

8c , 1 4 , 1 

m := -— (c 2 (1 + h + -) - 4), C := -(c 2 (l +h+-)-4), (so that b& + 2cC = 0); 
3b b 3 b 

2,1 2c 2 b , , 1 

^ t(c 2 (--2-2A)-1)(= — -6+-^l = 2(l-c 2 -c 2 / ; ) + -C). 
3 b b 2c 2 



Since 



1 1 wc/i / 1 \ ln*(c/i - f) _ ln*(-0 / mc/i / 1 



(t-ch) m t m t m+l \t m+2 )' (t-ch) m t m \ t \fln(-f) 
at t — -oo, we find that 

2c 2 /fc ytln(-f) 4c 3 / ; /fc rin 2 (-Q / i 1\ /ln(-f) 
0r(f-cA) := -^ + —^ + —^ + —^4 [ 1+0 ^n + 



Then a straightforward computation shows that 



Ri(t) := 0£(f) - c«^(f) + cf> T (t)(ifr Q (_t) - M0) = r n + r n 
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ln 2 (-f) ln(-f) 0(1) 



where 

r n := 2cT + 2c 2 Q/b + r n := 12.71 + 3\F - 4Jlc 2 /b - 2cT; 



and R 2 (t) := if/^t) - a/fat) + b<p T (t - ch)( \ - ^(f)) = r 2l + r 22 + 

with r 2 i := bT + 3cQ, r 22 := 6C(1 - c 2 ) + 3fec/i^[ - 2cg. 

Lemma 21. There exist T, Q and cr — cr(T, Q, c, b,h) > e such that <pr(f) > 0, tf/Q(t) > and 
R x (t) < 0, R 2 (t) < Ofor all t < -cr. 

Proof. Take T, Q such that r\\ > and r 2 \ < 0. Then it is easy to see that there is cr = 
cr(T, Q, c, b, h) > e such that t 2 <f> T (t) > 0, f 5 /?i(f) > and ttff Q (t) < 0, t 4 R 2 (t) < for f < -cr. □ 

Lemma 22. For every e > f/iere are r„, (?„ sufficiently large in absolute value and cr„ — 
o~(T n , Q„, c, b, h) < —e, such that the above defined functions 4>t„, >Pq„ 

I. are positive and strictly increasing on the interval (— oo, cr„]; 
//. are strictly decreasing on the interval (o~ n , cr„ + ch]; 
III. <f>' Ti (cr n ) = t[>' Qn (cr n ) = and 4> T Jcr„) + i// Qll (cr„) < e; 
TV. R\(t) > and R 2 (t) > Ofor all t < cr„. 

Proof. Take*-,, e [0.34,0.98] c (1/3, 1) and consider the sequences T„ — > -oo, Q n = -K n bT n jc — > 
+oo. It is easy to see that r\\ < and r 2 \ > for all sufficiently large n. 

Now, it is clear that, for a given fixed interval [z, —e], we have that <pr{t) < 0, t e [z, -e] for all 
sufficiently large negative T. On the other hand, for each T, function <pj is positive and strictly 
increasing on some interval (— oo, v]. These simple observations show that to every positive e 
we can indicate To < such that, for each T < To, the functions tpr, <P' T are positive on some 
maximal interval (-oo, cr\(T)) and <f>j(cri) < e/2, <f>' T (cri) = 0. The equation (p' T (cr\) = can be 
written as T — T\ (cr) with Pi satisfying 

C 2 cr] 

ri(cn) := -- 1 (1 + o(D), o-i -» -oo, (21) 
o ln 2 (-o"i) 

and strictly increasing on some maximal interval (— oo, y\(c, b, A)]. Hence, we see that crj = 
cr^r) depends continuously on T and monotonically converges to -oo as T — > -oo. 

Furthermore, since the equation r = F^o - !) has only one root cr\ e (-oo, y\(c, b, h)], we may 
suppose that (p' T {cr) < for all cr € (cr\ , 0"i + ch]. 

Using (fJTJ and the monotonicity of T\, one can readily establish that 



V=TMn(-r) 

o"i(r) = (i + o(i)), t -oo. 

2c 

Similarly, there is Qo > such that, for each Q > Qo, the functions \j/Q, tfr'g are positive on some 
maximal interval (-°°,cr 2 (Q)) and i//q(ct 2 ) < e/2, ifrQ(cr 2 ) = 0. Equation if/'gicr^ = can be 
written as Q = T 2 {cr 2 ) where 

r 2 ((T 2 ) = ~ ^ 2 - (1+0(1)), (T 2 -» -oo, 

3 In (-cr 2 ) 
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strictly decreases on some maximal interval (-oo,y 2 (c, b, h)]. From this we deduce that cr 2 = 
o~2(Q) depends continuously on Q and monotonically converges to -oo as Q — > +00. Also we 
may suppose that i[f' (cr) < on (o"2, <X2 + c/i]. Next, we have that 



criiQn) = - yj^QnQnXl + o(D) = o- x (T n ) yj ^(1 + o(l)), n ^ +», 

and since 3k„/2 e [0.51, 1.47], it is always possible to choose k„ in such a way that cr\{T n ) = 
ciiQn) '■- °~n for all large n. Obviously, k„ — > 2/3. 

Next, taking Q = r = 7 1 ,,, we find that for some functions a j, f3j, uniformly on n 
satisfying a } (t) = o(l),j8/(f) = (9(1), f -> -00, it holds 

ln 2 f— t) k bc~ l T 2 
Ri(t) = {2cT n {\ - K n + ai(t)) +A(0)— ^ - " 7 " ln 4 (-f)(l + «a(0). 

In this way, we prove the existence of 5] < -e which does not depend on n and such that R\{t) < 
for all t from some fixed interval (-00, 5\\. Thus we may assume in the sequel that cr„ < 6\ . 
Analogously, we can use the representation 

Ri(t) = (bT n (l - 3k„)) + a 3 ( f ))^-^ + KnbC f T " ln 4 (-f)(l + a 4 (t)), 

to establish that R2 is positive on some maximal interval (-00, 5 2 ), where 62 = 62(11) depends 
on n, lim^2(n) = -°° and /?2(#2(«)) = 0. Analyzing the latter equation, we find that there is a 
sequence b n — > b such that 



3c T n ln 2 (-6 2 (n)) f ... J2bT n 2/3 

— = -r— , and therefore 6 2 (n) = J —— \n l/i (-T n )(l + o(l)). 

2b„ 6 3 2 (n) V 27c 

Again, we have that <x„ < 62(11) for all large n, so that, without restricting the generality, we may 
suppose that both /?i(f),/? 2 (f) are positive on (-00, cr„], □ 

Remark 23. For r = 1 and small positive A, we will define <p A , <Pa by 

4> A (t) = Mt-A- 1 ), tMO = iPgit-A- 1 ), 

where sufficiently large T, Q are chosen as in Leinma^J] It is clear that for every To 6 R there 
exists A > such that cp A (t) > 0, faff) > 0, <p' A {t) > 0, ift' A (t) > Ofor all t < r , A e (0, A ]. /n 
addition, lim A ^o+((p A \t), ^ A (t)) = (0, 0), k = 0,1,2, uniformly on (-00, to]. 

Now we are in the position to prove Theorem [P7I 

Proof. Consider the functions 

+ (f,A) = min{l,0 A (O + cf> + (t)}, V + (t,A) = minjl,^) + *A+(0), 

where (p A , 4>a are defined in Remark |231 if r — 1 and in Lemma [T9l for r e (0, 1). By Remarks 
l20l 1231 and the implicit function theorem, there exist smooth functions 11(A), 12(A), such that 
lim A ^o+ ti(A) = fi,lim A ^ + 12(A) = t 2 and 

<D + (f,A),<D;(f,A) > 0, t < i 2 (A), with <D + (f,A) = 1, r > i 2 (A), 
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^+(t,A), x i'' + (t,A) > 0, t < ti(A), with W + (t, A) = 1, f >n(A), 

<D;( (2 (A)-,A) > <D' + ( t2 (A)+,A) = 0, r + ( tl (A)-,A) > <F' + ( tl (A)+,A) = 0. 

We claim that for t + i\(A), 12(A), d\, . . .cIm, and for sufficiently small positive A, the functions 
+ (f) := <D + (f,A), ¥ + (f) := <I> + (f,A) satisfy the system 

r 0;'(f) - c<D' + (f) + <D + (f)(l - r - <D + (t) + r¥ + (f)) < 0, 

I ?" + '(0 - c¥' + (f) + M> + (f-c/z)(l -¥+(*)) < 0, (22) 

[ > <d;(^+), v+(d r ) > v+idj+l 

Since differential inequalities (l22l hold trivially for all f > 1* := max{ii(A), 12(A)) (when <£+(?) = 
*F + (f) = 1), it suffices to prove ( l22l i for f e (-00, 1*), We will consider the following three cases. 
Case I . Let t < 1* := min{ii(A), 12(A)), then by Lemmas [T91I2T1 for all small A > 0, 

T;(t) - cY + (t) + b® + (t - ch)(l - < C(0 " cKW + b<f>At ~ ch)(\ - iff At)) 

-b(<p + (t - cWM + Mt - c*W+(0) < 0, (23) 

due to assumption D2 of Definition [T5l and the positivity of <f> + , iff a, <\>a, iff+. In a similar way (but 
this time using assumption Dl) we can evaluate F defined by 

r : = <D';(f)-co;(f) + o + (f)(i-r-o + (f) + ^ + (f)) 

= <t>'l(i) - C( p' + (t) + + (f)(l - r - (p + (t) + nfr+(t)) + <(0 - c<p' A (t) 

+ A (O(1 - r - A (O + r^ A (0) - 2^(f)0+(O + r0 A (O<MO + ^+(0^a(0- 

If r e (0, 1) then we obtain 

r = AO(e (J+v) ') + A 2 0(e' ,( * +1) ') + 0"(f) - c#.(f) + + (f)(l - r - + (f) + nft + (t)), 

and if r = 1 then 

(-t") m e n 

r = 0(^— ) + Rtf - a- 1 ) + <'(f) - cf + (0 + <M0(-<M0 + <A+(0)- 

In each of these two cases, for some small Ao, we obtainT < C(-f) m e v '0f(v, c)+o(l)), A € (0,Ao], 
t — > -00. Thus there exists r» < z'» such that T is negative for all t < t„ uniformly on A e (0, Ao]. 
On the other hand, since liiru^o min{ti(A), 12(A)) = min{fi,?2), we deduce from D2 and the 
above asymptotic representation of T the existence of A\ e (0, Ao) such that T is negative for all 
t G [t«,z'»], A € (0,Ai]. Thus (1221 holds for t e (-00, t*] and sufficiently small A e (0,Ai). 

Case II . Suppose now that i\(A) > l 2 (A) and let t e t*] = [12(A), 11(A)]. We have 

<D;'(f) - cO' + (r) + <D + (f)(l - r - <D + (f) + rT + (0) = -Ki - ¥+(*)) < 0, 

and, for sufficiently small A, ^'(f) - c^(f) + M> + (f - ch)(l - W + (t)) = 

Y := W1(t) - c^'Jt) + b((f> A (t - ch) + + (r - ch))(\ - ¥+(*)) < 0, f e [1., t. + eft]; 
l^ + '(0 - q^ + (f) + - + f A (i) - cif/' A (t) - bt// A (t) < 0, te[u+ ch, 1*]. 

Here we recall that i[i"(t) - ci//' + (t) +b(l— (M0) is negative on [?2 + ch, t\\ due to assumption D4. 
On the other hand, by the same assumption, we have that, for t € [u, 1* + ch] and all small A, 

T = ifr'Kt) - cifr'Jt) + bcf> + (t ~ ch){\ - ^ + (f)) - KM* ~ <*)^a(0 + 0a(* - <*)^ + (0) 

+^'(0 - af/jSf) + fe0 A (r - <*)(1 - < 0. 
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Case HI . Similarly, if l\(A) < 12(A) then for t e [«i(A), 62(A)], we obtain 

0>;'(f) - cO' + (0 + O + (0(l - r - <D + (f) + r¥ + (0) = ^(O - cO;(f) + <D + (f)(l - + (f)) 
= - C(j)' + + + (l - <p+) + <p' A ' - C(j)' A + <p A (l - (j) A ) - < 0, 

for all small A. Additionally, W^(t) - cW + (f) + b® + (t - h)(l - Y + (r)) = 0. Since <b' + (dj-,A) > 
<3)' + (dj+,A), x ¥' + (dj—,A) > y ¥' + (dj+,A) are obviously true for all small positive A, inequalities 
(l22l are proved for small A > 0. 

So let us fix some small A* > such that + (f) := <D+(r,A*), ¥+(0 := (!>+(?, A*) satisfy (l22l . 
In the continuation, we will prove the existence of lower solutions O : R — > [0, 1), which 
are defined as smooth non-decreasing functions satisfying the following system: 



We will treat separately each of the following cases: r e (0, 1) and r — 1. 
Suppose that r e (0, 1). It follows from the definition of <5 + (f, A*), < ¥+(t,A*) that for some 
positive ki = £i(A*), £ 2 = h(A*), 



Set now *?_(/) = 0, and define <&_(?), ? e R, as a unique (up to a translation) traveling front 
solution of the KPP-Fisher equation 



It is well known [l that <3>_(f) is strictly increasing and that <D_(f + so) = 0.5fci(A*)e'" + 0{e {A+5)t ), 
t — * — 00, for some small positive 5 and for an appropriate shift so = Sq(A*) which can be supposed 
to be zero. Hence, as a consequence of all mentioned properties of <1> ± , V I' ± , r e (0, 1), without 
restricting the generality, we may further assume that the third inequality in (l24l i is also satisfied. 

Let now r — 1 . For sufficiently large re (such that 0r„(o"n) < 1) | Ae„( cr n) < 1)> we consider the 
following C 1 -smooth increasing functions 



Take now re sufficiently large to have T„ < T, Q n > Q and <x„ < -(A*) l . Then 

4> A ,(t) = Mt - (A*)' 1 ) > <t>T„(t + O-n), Mt) = ip Q {t - (A*)' 1 ) > (f> Qn {t + cr„), t < 0. 

As a consequence, 




(24) 



<t> + (t,A*) > k x e At , T + (f,A*) > k 2 e M , t < 0. 



(p"(f) - ccp\t) + 0(f)(l ~r- 0(f)) = 0, 0(-oo) = 0, 0(+oo) = 1 - r > 0. 




Lemma [22] then implies that, for all f e R, 



0>'_'(f, re) - cQ>'_(t, n) + 0>_(f, «)(¥_(*, n) - 0_(r, re)) > 0, 
¥l'(f, re) - c*P'„(f, re) + fc«D_(f - ch, re)(l - ¥_(*, re)) > 0. 



<D + (f,A*) 
V + (f,A*) 



min{l,0 A .(f) +0+(f)} > 0r„(f + cr«) = 

min{l,0 A .(f) + (^+(0) > <Ar„( f + °~«) = T*_(f,re), t e R. 
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In order to finalize the proof of Theorem [TT] for a fixed negative number B < — (1 + r + b), 
we consider nonlinear operators 



tm, m) = 0(O(i -r-B- <p(t) + r<A(0), nit, m) = m - c/oa - m) - mt)- 

It is easy to check that T\,Ti are monotone in the sense that Tjdfii, ft\)(f) < Tj{<p2, 1A2HO, t 6 R, 
if < <f>i(t) < cf> 2 (0 < 1, < iAi(f) < ^(O < 1, t e R. Let zi < < zi be the real roots of 
the equation z 2 - cz + B — 0. Then every bounded solution (<p, iff) of differential equations in © 
should satisfy the system of integral equations 

<P(t) = N\ (0, m\ Mt) = N 2 {<t>, m\ where (25) 
Njty,ifr)(t) :=—!—( f e Zi( -'- s) Tj(cp^){s)ds + f e z *- s) Tj(<f>,ft)(s)ds 

Zl - Zi \J-00 Jt 

Conversely, each positive strictly monotone bounded solution (<p, iff) of $25[ yields a wavefront 
for {8). It is clear that the operators Nj are also monotone. Additionally, it is easy to see that 
Nj{(f>, tAXO is increasing if both <f>, if/ : R — > [0, 1] are increasing functions. 
Hence, taking into account (1221) . d24l) and Lemma[l8] we conclude that 

*-(*) < O a) (0:=M(1>-, v l'-)(f)<M( ( l> + , ( I' + )(0:=O ( H 1) (r)<O + (f), 
<F_(?) < ^(f) := AfcC*-, <F_)(f) < AT 2 ($ + , ¥+)(*) := ¥^(0 < T\(f). 

Therefore the sequences of positive uniformly bounded (by from below and by 1 from above) 
monotone continuous functions 

>pC ,+1 >( f ) = Ar 2 (®?\ V^XO. 4>? +1) (0 = M(0 W , ^XO, n = 1,2, . . ., (26) 
and ¥? +1) (0 = A^O^^Xf), *? +1) (0 = MC^.^KO, «= 1,2,..., 

are strictly increasing and decreasing, respectively. Set <I> = lim®^', = lim v I' < " ) , then 

(D < O < <D + , >T < ¥ < T+. (27) 

Furthermore, a direct application of the Lebesgue's dominated convergence theorem to (l26l 
shows that the pair (<I>, *P) solves system (1251 ). Since <t>(f) > for all t, we may conclude from (l25T l 
that ^(f) > 0, ? e R. Note also that <5(-oo) = ¥(-00) = in virtue of (J27]). Now, since <5, *F are 
positive, increasing and bounded functions, the values of 0(+oo), v I'(+oo) are finite and positive. 
A standard argument based on the Barbalat lemma (cf. [27]) shows that <t>(+oo) = l F(+oo) = 1. 

Finally, the validity of asymptotic formula (0 follows from d27l i. To prove ©, we first 
observe that, due to ( l27T i and Theorem|6] there exists T/, < such that O(f) and ^(f) are bounded 
(from below and from above) by cie M for some c, > and all t < T/,, Then we can apply 
Proposition 7.2 from 111 811 to the first equation of (|8} in order to obtain the desired formula for 
<!>(?)• Using this formula and the change of variables y = ft - be x{, ~ ch) j '(1 - r), we then get easily 
the second formula of ©, cf. the proof of Lemma[l4]and Corollary[T2] □ 

5.2. Proof of Theorem\7\ 

The simplest form of regular super-solutions <f> + , ft + is exponential, we can write them as 

4>+(t) = e v \ ft + (t) = De vl , D = be- vch /(cv - v 2 ), (28) 

where D is chosen in such a way that the second inequality in D2 as well as D4 were satisfied 
(details are given below). In order to simplify the notation, in the sequel we will write b' = be~ vch . 
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Lemma 24. Suppose that v + jA is close to c/2, r > 0, and 

(cv- v 2 )(l + 1/b') > 1, cv- v 2 + b'. (29) 
Then ( I2<SD determines a regular super-solution for (0. 

Proof. Clearly, Dl is satisfied with C\ — 1, t\ — v _I ln[(cv- v 2 )/b'] + 0, f 2 = 0. Still we have to 
check hypotheses D2, D3, D4. Depending on the sign of t\, we will analyze the next two cases: 
Case I. t\ > = ? 2 or, equivalently, < b'/(cv - v 2 ) < 1. If t e [0,*i], then D4 holds 
because 

f ~b'e vt i// + (f) <0, f e [0,c/i]; 

^ + (O-c^ + (O + 6fflin{1.0 + (*-^)Kl-MO) = ^yW) < 0, r e [cMil 

If f < 0, we have that 

if/lit) - cip' + (t) + b<f> + (t - ch)(\ - <p+(t)) = -b'De 2v ' < 0; 

- ccf>' + (t) + + (O(l - r - + (r) + /^+(0) = e v 'Mv, c) - + (O + r^ + (t)J < 0. (30) 

Case II. Next, let f] < f 2 = so that b' [(cv - v 2 ) > 1. If t e [ft, 0] then 

<(f) - c<(0 + + (O(l - + (O) = e v '(l - e v < + v 2 - cv) 
< e v '(l - e vh + v 2 - cv) = e v '(l - (cv - v 2 )(l + 1/b')) < 0, 

and D3 holds. Now, for t <t\, condition D2 is true since 

#(0 -c#(t) + 0+0(1 -r-* + (0 + ^+(O) = (31) 
e v '(x(v,c) + e v '(-l + rb'Kcv - v 2 )) < e v 'maxtr(v,c), 1 - (cv - v 2 )(l + l/b')} < 0; 

^ + '0 - c^(0 + b<f> + (t - ch)(l - ifr + (t)) = -b'De 2vl < 0. 

This completes the proof of the lemma. □ 

Corollary 25. The existence statement of Theorem%7\holds true. 

Proof. First, we assume that c > c#. Then clearly there is a positive v meeting all requirements 
of Lemma l24l This assures the existence of a regular super-solution for (0. By Theorem [P71 
system ^ has a positive monotone wavefront. 

Next, we consider the case when c = c#, r € (0, 1], Let cj I c# be a strictly decreasing 
sequence of velocities and (<f>j, tj/j) be a sequence of corresponding traveling fronts (existing in 
virtue of the first part of the proof). Since 

= (pj(-oa) + iffj(-oa) < (f)j(t) + iffj(t) < (pj(+oa) + lffj(+oo) = 2 

and the function <f>j(t) + i/fj(t) is increasing in t for each fixed j, we may assume that <f>j(0) + 
i/fj(0) = 3/2, j — 1, 2, 3, . . . Using the standard compactness arguments and then applying the 
Lebesgue's dominated convergence theorem to the system of integral equations (l25T i: 

4>j(t) = Mi(<pj,i//j,Cj)(t), ifr/t) = N 2 (4>j,ifrpCj)(t), 

we may assume, without restricting the generality, that lim j(<f>j, iffj) = (0, 0-) uniformly on bounded 
intervals, where (0, 0-) is a monotone solution of ([8} with c = c#. Since (0, 0)(±oo) are steady 
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state solutions of (O and 0(-oo) + 0(-oo) < 0(0) + 0(0) = 3/2 < 0(+oo) + 0(+oo), we find that 
necessarily 

0(-°°) = 0, 0(-°°) e [0, 1], 0(+oo) = 0(+oo) = 1, 

(if 0(-oo) > 0, then 0(-oo) = 0(-oo) = 1 and thus 0(0) + 0(0) = 2, a contradiction). To finish 
the proof of the corollary, we have to establish that 0(-°o) = 0. In order to prove this, we can 
apply the part [A] (for r e (0, 1)) and the part [C] (when r — 1) of Theorem|6]to find that either 
iffj(i) < K(f>j(t), t e R (for r e (0, 1)) or y (f) < ^Mfi/t), t e R (for r = 1). Therefore either 

0(0 < K<t>(t), t e R, or 0(f) < y/M<f>(t), t e R, so that 0(-°°) = 0. □ 

5.5. Proof of Theorem^ 

Let now c < c# so that for some v close to c/2 

< (cv - v 2 )(l + — ) < 1, cv - v 2 * v * jA. (32) 

Analyzing the proof of Lemma [24] under these assumptions, we see that t\ < = ?2 and the 
main obstacle to develop successfully the proof of Case II appears when we want to estimate 
expression (IBTb near t\ < 0. Therefore we may expect that, after an appropriate modification 
of super-solutions (l28T l in some neighborhood of t\, the result of Theorem Q can be improved. 
Below, we develop this idea by considering 0+(f) = e vl , t e R, and C 1 - smooth function 

Mt)={ DeVt ' Vt-'r (33) 

\ p + qt, if f > f„. v y 

Here p, q, U will be chosen to satisfy the first inequality in D2 for all ( e R. The mentioned 
inequality can be written as 

+ (f) < y(f) := r~\cv- v 2 + r- 1 + e vf ). (34) 

Now, assuming ( 1321 ) and analyzing the mutual positions of convex graphs of the functions y(t) 
and De VI , we deduce that these graphs should have exactly one point of intersection (or tangency) 
below the level y = 1 . Indeed, otherwise cv-v 2 + r-l > implies that y(t) > De vt for all t where 
y{t) < 1. As a consequence, 1 = j{sq) > De vs " at some sq which implies d29l ), a contradiction. 

The above consideration and a direct computation show that there exists a a unique line 
y — p + qt which is tangent to the graphs of De vt and y{t) at the respective points ?* < f . From 
the tangency conditions q = Dve vt " - y'(f ), p = De v '" - qt„ = y{f ) - qt*, it follows easily that 

q Dve (cv - v 2 + r - l)v 1 q 

p = — In , q = , t t — — In — . 

v q r\n{rD) v Dv 

It follows from the above construction that + defined by ( l33l is C 1 - smooth and 0+(f) < y(t) for 
all t + f . It is clear that, after making an arbitrarily small change of p, q, t*, we may assume that 
<A+(0 < r(0 for all t 6 R. 

Hence, taking + as in (1331 and 0+(f) = e vr , we have to check only the second inequality in 
D2 on the interval [f t , +oo). This inequality can be written as b'(l - p)/q < b't + ce~ vl . Since 
y = b't + ce~ vt has a unique critical point (an absolute minimum) at f = v~' ln(cv/Z?'), the latter 
inequality amounts to 

(1 - p)v < q]n(ecv/b'). (35) 
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After recalling the definition of p and D and taking into account that v can be chosen as close to 
c/2 as we want, we rewrite (l35l l asw< 2(2 + lnw), where 



2r Ab'r c 

c z /4 + r - 1 c 1 q 

Notice here that the assumptions c 2 > 4(1 - r) and c 2 < 4/(1 + (Z>') _1 ) imply r(fo' + 1) > 1 and 
c 2 < 4b' r so that w > 0. Furthermore, since a> is decreasing in c 2 /4, we find that 

2r , 2r(\+b') 

oj > In rib' + 1) = — In rib' + 1) > 2. 

b'/il + b') + r- 1 r(l +/?')-! 

A direct graphical analysis shows that the interval u> e (0.14555 . . . , 8.21093 . . .) gives the solu- 
tion of oj < 2(2 + In cj). In consequence, since we additionally have u > 2, the latter inequality 
is equivalent to to < io# = 8.21 . . . which can be written as (0. This proves Theorem|8] □ 



6. Proof of Theorem|9] 

The proof is divided into three claims. 

Claim I : The propagation speed c* is unique. Indeed, suppose that (0i, tff\, cj) and i<p2, i^ 2 , ci), 
c\ < ci, solves the nonlinear eigenvalue problem ©. It follows from Lemmas [T3l and [T4l that 
there exist p < q such that ii\it) > iffiit) for all t e R \ [p, g]. As a consequence, the closed set 

S := {s : (Ai(f + j) > 1A2W, teil^K 

is non-empty and has a finite s* := inf»S. It is clear that i//\it + s*) > fait), f e R, and since 
always ij/\it + s*) > t^ 2 (f) for t <k -1 and t » 1, we deduce that (t + s«) = ^(t) at some point 
r (otherwise > inf >S). By a similar argument, there exists t, such that 0i(f + f») > (fait), t e R, 
and 0i(T + f«) = </> 2 (r) for some T. Suppose first that t m < s m . Without restricting the generality, 
we may assume that s* = 0, r = 0. Then £, < so that </>i(f) > <piit), t e R, and thus we get 

= OAi - <A2)"(0) - cityi - i/f 2 )'(0) + ic 2 - ci)^ 2 (0) + bifai-dh) - $ 2 i-c 2 h))i\ - (AK0)) > 0, 

a contradiction. Next, suppose that t t > s t . We may assume again that that f* = 0, T — 0. Then 
5* < so that i^i(f) > il/lit), f G R, and thus we get 

= (0i - 2 )"(O) - ci(0, - 2 )'(O) + (c 2 - ci)^(0) + r0!(O)OAi(O) - <A 2 (0)) > 0, 

a contradiction. Hence ci = c 2 and Claim I is proved. 

Claim II : c* < c m :- min{c#, c a }. Let (<p t , ifr*, c+) be the solution of ([8j. On the contrary, suppose 
that c* > c,„ and take an arbitrary c' e (c m , c*). Then (0„, c') is a lower solution: 

(p"it)-c'<f>',it)+(f> m it)il-r-(p m (t)+ri// m it)) > 0, iy,'if)~cY m if)+b(p m it-c'h)il-iJr,it)) > 0, r e R. 

For the same c' we consider the upper solutions <D+(f) = min{l, </>+(f)L = min{l,i/r + (f)}, 

with <f>+,i//+ defined in Subsections 15.21 IB31 By Lemma [T3l we may suppose (possibly, after a 
translation of (0„ tf/ t )) that 0,(f) < 3>+(r), <A*(0 < v P+(0> f 6 K - But then there exists (cf. the last 
part of Subsection l5.ll starting from formula (1231 ) a monotone traveling front propagating at the 
velocity c' < c*. However, this contradicts to Claim I. 
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Claim HI : Set c+(/i) :— c*(r, b, h)for some fixed r,b > 0. Then c+(h) is a non-increasing function 
on its domain. Suppose that c+ihx) > c+ihf) for some h\ > ho. Let (<pj, c+ihj)) be respective 
solutions of (|S). Then, for a fixed c e (c+ihf), c*(/zi)), it holds 

0' 1 '(O-c0i(O+0i(O(l-> , -0i(O+r^i(O) > 0, ifriW-af/^+bfrit-chiXl-ij/iii)) > 0, t s R, 

^(t)-c4>' 2 (t)+tfr 2 (t){l-r-4>2(t)+rifr 2 (t)) < 0, ^'(f)-c^(0+^ 2 (f-c/n)(l -^(0) < 0, f e R. 

Moreover, due to Lemmas [T3l and fT4l we may assume that 0i(f) < fait), < 412(f), (el. 

Therefore (<pj, if/j, c), j - 1,2, forms a pair of upper and lower solutions for (|8} considered with 
c and h\. As a consequence, system ((SJ with /1 = hi has two different propagation speeds: c and 
c*(/zi) > c. However, this is a contradiction with Claim I. □ 
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